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SHORT COMMUNICATION 

The utility of ITS2 in spider phylogenetics: notes on prior work and an example from Anelosimus 
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Abstract. The ribosomal internal transcribed spacer ITS2 is probably the most popular nuclear DNA marker used to 
examine relationships among and within species in animals and plants. ITS2 sequences have also begun to be used as DNA 
barcodes. ITS2, however, has rarely been used in studies of spiders. Here, I examine the potential utility of this marker for 
spider phylogenetics based on preliminary data for Anelosimus spiders and a brief summary of prior work. The secondary 
structure of ITS2 facilitated alignment of highly divergent sequences and indicated that secondary structure morphology 
might be phylogenetically informative in itself. Phylogenetic analysis of Anelosimus species was congruent with a prior 
study based on a combination of six mitochondrial and nuclear loci plus morphology regarding the deeper clades within the 
genus. However, ITS2 had insufficient variation to resolve relationships within species and among closely related species. 

Previous studies have also discovered relatively little within-species variation in ITS2. In sum, ITS2 is an easily amplified 
and sequenced marker that is underutilized in spider phylogenetics; however, it has limited uses at the lowest taxonomic 
levels and is not likely to be a universally useful DNA barcode marker. 

Keywords: Phylogeny, DNA barcode, ITS2 secondary structure, sociality, Theridiidae 


ITS2, which is flanked by the 5.8S and large ribosomal subunit 
(28S) nuclear genes, is perhaps the most popular marker used to 
resolve relationships among*and within species in animals and plants 
(Alverez & Wendel 2003; Bailey et al. 2003; Young & Coleman 2004; 
Schultz et al. 2005; Coleman 2009; Schultz & Wolf 2009). ITS2 
sequences have also been proposed as effective DNA barcodes (e.g., 
Ben-David et al. 2007; Park et al. 2007). The popularity of this marker 
stems from a generally high level of variation, yet relatively conserved 
secondary structure, and ease of amplification and sequencing. 
However, comparatively few studies on spider phylogenetics have 
utilized this marker despite these benefits and a general paucity of 
good primers for nuclear markers. Among the few ITS2 studies in 
spiders, most focus on low taxonomic levels, reconstructing relation¬ 
ships among, and in some cases, within, species (Hedin 1997; 
Hormiga et al. 2003; Arnedo & Gillespie 2006; Chang et al. 2007; 
Bond & Stockman 2008). In spiders, ITS2 has generally been found to 
be a useful marker offering resolution at the species level, especially so 
in more genetically structured systems such as in trapdoor spiders 
(Bond & Stockman 2008), cave dwelling nesticids (Hedin 1997), and 
island radiations (Hormiga et al. 2003; Arnedo & Gillespie 2006). 
Other studies have used ITS2 as a tool to help separate closely related 
species. Variation allowing separation of closely related species/ 
populations was found in Poltys (Smith 2006), Parclosa (Chang et al. 
2007) and Latrodectus (Vink et al. 2008). However, variation was 
insufficient to separate closely related North American Latrodectus 
species (Zhang et al. 2004) or populations of L. katipo Powell 1870 
(Vink et al. 2008). 

This note reports on the utility of ITS2 data to resolve phylogenetic 
relationships among and within Anelosimus spider species, well 
known for their multiple origin of social behavior (Aviles 1997; 
Agnarsson 2006; Agnarsson et al. 2006). I use exemplar species from 
across the phylogeny of the genus and specimens representing most of 
the known 16S mitochondrial haplotype diversity within one species, 
A. eximius (Keyserling 1884). As the goal of this paper is practical 
application, I do not see a reason to prune the analyzed matrix to the 
exact ITS2 sequences, but I refer rather loosely to the entire region 
amplified by the FITS and RITS primers (see below) as ITS2. 

I collected specimens in the field and placed them in 95% ethanol. 
Genitalia were abscised and stored as vouchers at the Zoological 


Museum of the University of Puerto Rico, while DNA was isolated 
from each individual using the prosoma, the abdomen, or both, with 
the QIAGEN DNAeasy Tissue Kit (Qiagen, Inc., Valencia, CA). 
I used the ITS-5.8S (FITS) and ITS-28S (RITS) primers (White 
et al. 1990) (FITS GGGACGATGAAGAACGGAGC, RITS 
TCCTCCGCTTATTGATATGC), using standard protocols with 
an annealing temperature of 47° C for 30 cycles. The PCR products 
were sequenced by the MACROGEN service, and sequences were 
submitted to GenBank (Accession numbers: HM584843-HM584883). 
Data for outgroups ( Latrodectus , Enoplognatha), were obtained from 
GenBank. Preliminary alignments were done using ClustalW 
(Thompson et al. 1994) with gap opening and extension costs set at 
24/6 and 8/2. Most of the sequences aligned readily; however, these 
preliminary alignments revealed an area of a particularly difficult 
alignment. Analyses of Clustal aligned matrices gave results largely 
incongruent with prior phylogenetic hypotheses, which were based on 
more data, mostly due to the placement of the root of the Anelosimus 
tree. Therefore, preliminary alignments were followed by manual 
and automated alignments taking into consideration the implied 
secondary structure of ITS2 (Fig. 1). I used the ITS2 database (on¬ 
line at http://its2.bioapps.biozentrum.uni-wuerzburg.de/cgi-bin/index. 
pl?about) to annotate the sequences and find the 5.8 and 28S flanking 
regions. Anelosimus ITS2 sequences were generally short, ranging 
from 223-305 bp. Non-ITS2 sequences were then removed and ITS2 
secondary structure implied using the 4Sale software (Seibel et al. 
2006, 2008). A standard model of ITS2 secondary structure was 
developed by Sehultz et al. (2005). This model has one long arm, or 
‘helix’ (helix III), and three shorter helices (helices I, II, and IV), all 
four radiating from an area of a large loop. Secondary structure 
analyses reveal that the region that aligns poorly using Clustal 
corresponds to helix IV of the consensus ITS2 secondary structure 
model (Schultz et al. 2005) which is present in Latrodectus , A. 
rupummi Levi 1956, and a short version of it in Enoplognatha , A. 
nigrescens (Keyserling 1884) and A. ethicus (Keyserling 1884), but lost 
in the ‘ eximius lineage’ (Fig. 1). 

Once this helix is identified, manual alignment of this region is 
facilitated, essentially aligning apparently homologous regions of 
helix IV in those taxa that have it, and inserting a gap for the entire 
arm region in those taxa that lack it. Automated alignment was also 


377 


378 


THE JOURNAL OF ARACHNOLOGY 


-A 

c c 


Co 

K ) 


1 C 

* * 


' c « 


«?-<? 

c«fc 

* t 


<? o 

*< 


O# 

\ g 

*eV 6 

V 


<?■ * 


£-<? 

<?<? 


K C ® I 

«o!o 1 




#*L 


O 


V*. 




g H$ c 6 *C C# <5 C* 


a g k c 

,t>\» 
*- § % 


i»\ »A C 

.5. ?*f. c 

>' 


<5M 


<r C A 

:J 


S } 
♦ ** 


I> 

«?#* 




* * <?<?•<? 

5 oc 

■4 f ® 

*'£& c 4 * 
oc 

* dC # 

*»*• 

<*<? 


XA <* 


«, c V* 

&v 


»*» 

G C 
G'% 

C G 
<16 
i>i 
6 # 
6G 


« V _ 
**•>,. 


g4S G 

•Og 


$ * 


a ► 

OG 


O* »C 6 

t it 

Sc^CG 

* 


III 




IV 


M 


A. eximius 


C J 

C C 

ce 
'# 
c & 

4& 

<t * 

A. baeza 4 * • 


#GC« 
G C 

*9 

G C 
>• 
:i 

C-G 

• * 
g a 


A. rupununi 


Figure 1.—Secondary structure (Bruccoleri layout) of ITS2 as implied by 4SALE for two species of the eximius lineage (A. eximius and A. 
baeza) and A. rupununi. The overall similar secondary structure reflects sequence similarity across most of the ITS2 sequence in these taxa that is 
readily alignable. However, A. rupununi has a helix (arrow), corresponding to helix IV of the ITS2 consensus structure of Schultz et al. (2005) that 
has been lost in the eximius lineage. A second region of difficult alignment is a loop region preceding this helix IV (stars). 


conducted with the 4Sale software, using the remote 4Sale option. 
The automated alignment was not modified other than by fixing the 
first eight aligned characters, representing a five base pair sequence 
identical in all taxa, which had been rather randomly spread out. The 
remainder of the automated alignment did not contain conspicuous 
areas of misalignment. Aligned matrices and results are available 
from the author upon request. 

The appropriate substitution model was selected with Modeltest 
(Posada and Crandall 1998), using the AIC criterion (Posada and 
Buckley 2004) with a parsimony tree chosen as the basis for 
Modeltest. The best model was GTR + F + I (Yang 1994). Bayesian 
analysis was performed using MrBayes V3.1.2 (Huelsenbeck and 
Ronquist 2001). The Markov chain Monte Carlo was run with four 
chains for 10,000,000 generations (repeated twice), sampling the 
Markov chain every 1000 generations, and the sample points of the 
first 5,000,000 generations were discarded as “burnin”. Maximum 
likelihood analyses were conducted in the program Garli (Zwickl 
2006), using the GTR + T + I model and 200 search replications. 
Parsimony analysis was done using TNT default settings under 
traditional search, with 1000 search replications. To calculate 
divergences among and within species in previous studies (Table 1), 
I downloaded sequences from GenBank via Mesquite (Maddison and 
Maddison 2009) and calculated uncorrected genetic distances in 
Mesquite. 

The phytogenies are largely congruent using the Bayesian, 
likelihood, or parsimony criteria, and whether based on the manual 
or automated alignment (Fig. 2); hence, only the Bayesian results are 
discussed. To the extent that the current results are comparable to 
prior studies that included more taxa, they recapitulate the deeper- 


level phylogeny of Agnarsson et al. (2007, 2010) based on six 
molecular loci combined with morphology (Fig. 2). However, the 
analysis does not resolve relationships among closely related species 
of the stucliosusijucuuclus groups and does not reflect strong 
mitochondrial population structure within A. eximius. 

Within Anebsimus , therefore, the utility of ITS2 seems very limited 
at the lowest taxonomic level (within species, between closely related 
species), but higher at intermediate taxonomic levels. Other studies of 
closely related theridiid species have also found little to no 
informative variation among closely related species (Zhang et al 
2004, Vink et al. 2008). However, in cases where population 
structuring is particularly strong, such as in trapdoor spiders (Bond 
& Stockman 2008) and cave-dwelling nesticids (Hedin 1997), ITS2 
was found to be useful at the interspecific, and even intraspecific level. 
Based on this and prior studies, ITS2 is a useful and readily 
obtainable marker for phylogenetic studies that look at relationships 
within genera and families of spiders. In general, intraspecific 
variation is low in spiders (about 1% on average, Table 1), as is the 
variation between sister species, but the variation differs across 
groups and is notably high in some trapdoor and cave dwelling 
spiders (Table 1). Closely related Anebsimus and Latroclectus species 
have very low interspecific variation (typically < 1%, about 0.7% in 
A. eximius , which shows high mitochondrial variation), insufficient to 
resolve relationships among closely related species, or to diagnose 
species. The utility of ITS2 at lower taxonomic levels thus will vary 
depending on the group. At higher taxonomic levels the main 
difficulty will be extreme sequence divergence (e.g., 27% between A. 
rupununi and A. eximius ), thus complicating alignment. However, 
ITS2 secondary structure can facilitate alignment of divergent 
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Latrodectus katipo 
Enoplognatha ovata 
A. rupununi 016A 
A. ethicus 021A 
A. nigrescens 019A 
A. dubiosusOMA 
A. jabaquara 062A 
A. domingo 009A 
A. domingo 109A 
A. eximius 132A 
A. eximius 133A 
A. eximius 135A 
A. eximius 202A 
A. eximius 138A 
A. eximius 127A 
A. eximius 203A 
A. eximius 233A 
A. eximius 232A 
A. eximius 126A 
A. octavius 026A 
A. octavius 015A 
A. octavius 027A 
A. baeza 028A 

A. oritoyacu 008A 
A. analyticus 023A 
A. baeza 018A 
A. fosu/r? 013A 
A. fosum 014A 
A. tosum 003A 
A. fosm7?024A 
A. baeza 001A 
A. guacamayos 002A 
A. guacamayos 010A 
A. guacamayos 007A 
A. sp. n. [viera] 020A 
A. studiosus 031A 
A. sfutf/'osus 011A 
A. studiosus 012A 
A. baeza 022A 
A. baeza 025A 
A. baeza 004A 
A. baeza 005A 
A. baeza 006A 


Figure 2.—50% Majority consensus from the Bayesian analysis of secondary-structure-informed manual alignment with numbers showing 
posterior probability support values. Major deeper level clades and species groups well supported by prior work are recovered: see clade labels. 
All labeled clades were recovered in all analyses, except the 6 domingo group’. However, at lower taxonomic levels very little variation was 
observed, resulting in low resolution. Within the ‘sclerotized CD clade’ (A. analyticus plus the jncundns!studiosus complex) only A. guacamayos 
Agnarsson 2006 and A. studiosus (Hentz 1850) were recovered as monophyletic, the relationships among species were largely unresolved and 
inconsistent with prior work (Agnarsson 2006, 2010; Agnarsson et al. 2007). Within A. eximius , a species showing population division and strong 
mitochondrial structuring, no phylogenetic structure was recovered. 














































380 


THE JOURNAL OF ARACHNOLOGY 


Table 1.—ITS2 maximum intraspecific sequence divergences, and estimation of maximium and minimum divergences between sister species, 
in previously published studies of spiders. Estimated intraspecific sequence divergence is likely conservative overall, as some species were sampled 
only by 2-3 individuals. However, even for species sampled by 10 or more individuals and from geographically distant localities (e.g., Latrodectus 
katipo , Anelosimus ex twins) the divergences were low. 


Family 

Genus 

Species or putative species 

Maximum intraspecific 
sequence divergence 

Reference 

Araneidae 

Poltys 

illepiclus 

0 

Smith 2006 

Araneidae 

Poltys 

stygius 

0 

Smith 2006 

Araneidae 

Poltys 

laciniosus 

0 

Smith 2006 

Cyrtaucheniidae 

Aptostichus 

clade 5 

0.039 

Bond and Stockman 2008 4 

Cyrtaucheniidae 

Aptosticlins 

Clade 2 

0 

Bond and Stockman 2008 

Cyrtaucheniidae 

Aptostichus 

Clade 3 

0.025 

Bond and Stockman 2008 

Cyrtaucheniidae 

Aptostichus 

Clade 1 

0.004 

Bond and Stockman 2008 

Linyphiidae 

Orsomvelles 

graphica 

0.005 

Hormiga et al. 2003 

Linyphiidae 

Or son we lies 

macbeth 

0.01 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

falstaffius 

0 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

polites 

0 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

othello 

0.005 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

amber sono rum 

0 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

malus 

0.002 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

arc anus 

0.02 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

calx 

0.002 

Hormiga et al. 2003 

Linyphiidae 

Orsomvelles 

vent ns 

0.02 

Hormiga et al. 2003 

Lycosidae 

Pardosa 

as tiger a 

0.03 

Chang et al. 2007 

Lycosidae 

Pardosa 

astigera (phenotype A) 

0.003 

Chang et al. 2007 

Lycosidae 

Pardosa 

astigera (phenotype B) 

0.005 

Chang et al. 2007 

Nesticidae 

Nesticus 

barri 

0.0025 

Hedin 1997 1 

Nesticidae 

Nesticus 

bar row si 

0.0102 

Hedin 1997 

Nesticidae 

Nesticus 

bishopi 

0.0051 

Hedin 1997 

Nesticidae 

Nesticus 

cooperi 

0.0059 

Hedin 1997 

Nesticidae 

Nesticus 

“ dellingeri 7 

0.0076 

Hedin 1997 

Nesticidae 

Nesticus 

gertsclii 

0.0152 

Hedin 1997 

Nesticidae 

Nesticus 

mimus 

0.0119 

Hedin 1997 

Nesticidae 

Nesticus 

nasicus 

0.0077 

Hedin 1997 

Nesticidae 

Nesticus 

silvanus 

0.0034 

Hedin 1997 

Nesticidae 

Nesticus 

stupkai 

0.0102 

Hedin 1997 

Nesticidae 

Nesticus 

carterr 

0.0321 

Hedin 1997 

Nesticidae 

Nesticus 

nov. sp 

0.0051 

Hedin 1997 

Nesticidae 

Nesticus 

paynei 

0.0076 

Hedin 1997 

Nesticidae 

Nesticus 

tennesseensis 

0.0085 

Hedin 1997 

Salticidae 

Havaika 

OK9, OK24, OW28, OW29 

0.0176 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

OK8, OK23, OW11L OW158 

0.03 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

H83, HI09, HI37 

0.005 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

EM 128, WM88, WM159 

0.018 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

K85, K86, K87 

0.03 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

EM81, MK82, WM89 

0.0025 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

H10, HI 10, EM90 

0.015 

Arnedo and Gillespie 2006 3 

Salticidae 

Havaika 

WM88, WM159, EM 128 

0.018 

Arnedo and Gillespie 2006 3 

Theridiidac 

Latrodectus 

katipo 5 

0.002 

Vink et al. 2008 

Theridiidae 

Latrodectus 

hasselti 

0 

Vink et al. 2008 

Theridiidae 

Latrodectus 

hasselti 

0.0027 

Zhang et al. 2004 

Theridiidae 

Latrodectus 

mactans 6 

0.014 

Zhang et al. 2004 

Theridiidae 

Anelosimus 

eximius 

0.007 

This study 

Theridiidae 

Anelosimus 

don lingo 

0 

This study 

Theridiidae 

Anelosimus 

tosum 

0.008 

This study 

Theridiidae 

Anelosimus 

stndiosus 

0.01 

This study 

Theridiidae 

Anelosimus 

guacamayos 

0.002 

This study 

Theridiidae 

Anelosimus 

octavins 

0.007 

This study 

Theridiidae 

Anelosimus 

baeza 

0.02 

This study 

Average 



0.01 
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Table 1.—Continued. 


Interspecific sequence divergence (sister species) 


Min 

Max 

Reference 

Nesticus 

0.40% 

~9% 

Hedin 1997 

Latrodectus 

0% 

0.50% 

Vink et al. 2008 

Latrodectus 

0% 

0.83% 

Zhang et al. 2004 

Havaika 

2% 

4% 

Arnedo and Gillespie 2006 3 

Pardosa 

2.50% 

6.70% 

Chang et al. 2007 

Orsomvelles 

0.70% 

5.90% 

Hormiga et al. 2007 

Anelosimus 

0.60% 

2.80% 

This study 

Poltys 

0.70% 

-10% 

Smith 2006 

Aptostichus 

2.20% 

5.30% 

Bond and Stockman 2008 


1 Note that multiple individuals within populations always had zero sequence divergence, interspecific sequence divergences reflect those among 
isolated populations 

2 Represented two species, each with intraspecific divergence less than 1.5% 

3 Informal species, reflecting putative species from Fig. 5 in Arnedo and Gillespie (2006), codes in ‘species’ column refer to specimens 

4 Sequences from ‘clade 4’ were not found on Genbank 

5 One variable site 

6 More divergence found within than between individuals 


sequences (Young & Coleman 2004) (Fig. 1). Based on my findings 

and those of Vink et al. (2008) and Zhang et al. (2004), ITS2 does not 

emerge as a suitable ehoice of universal DNA barcode. 
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